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Spin decay of ‘ Aryabhata ’ 


S K SHRTVASTAVA and C K RAJASINGH* 

Department of Aeronautical Engineering, Indian Institute of Science, 

Bangalore 560 012 

•Mission Operations and Planning Division, ISRO Satellite Centre, Bangalore 560 058 
MS received 24 August 1977; revised 7 August 1978 

Abstract. The spin decay of Aryabhata due to eddy current losses in the earth’s 
magnetic field is analysed in this paper. The satellite is idealised as a conducting 
spherical shell while the earth’s field is taken to be equivalent to that of a magnetic 
dipole. The results of this simplified analysis match quite well with the actual 
observations. 
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1. Introduction 

Aryabhata, the first Indian satellite, is spin-stabilised for inertial orientation. The spin- 
rate is maintained within a range (designed 90-15 rev/min, revised 52-5 rev/min) by 
occasional firing of cold-gas jets (Goel et al 1978). The mission life beinglimited essen¬ 
tially by the available gas aboard the satellite, it is important to analyse the spin decay. 

Some of the possible sources of spin-decay (or rise) are eddy currents, residual mag¬ 
netic moment, skin friction in rarefied atmosphere, structural flexibility, leakages, 
hysteresis, changes in mass and moment of inertia and solar radiation pressure. 
For Aryabhata structural flexibility can be ignored; magnetic cleanliness is ensured to 
a high degree; changes in mass and moment of inertia, due to gas depletion and defor¬ 
mations and movements in the tape recorder are quite small; the orbit is sufficiently 
high (about 600 km near-circular) for the atmosphere to affect the satellite; and the 
solar radiation pressure torque is negligible. However, eddy current losses due to 
the aluminium alloy structural shell spinning in the earth’s magnetic field cannot be 
ignored. A comparison of the predicted decay due to eddy currents (Rajasingh 1974 
unpublished) with the actual observations confirms this view. This paper presents 
the analysis and comparison of the results with the observations. 


2. Analysis 


Torque due to eddy-currents on a spacecraft rotating in the earth’s magnetic field 
is given by 


T = /"j» x (JxB )dV, 

•A list of symbols appears at the end of the paper 


( 1 ) 
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where J — a {v® +i(^xB)x ~py. (2) 

Potential $ is required to satisfy the boundary conditions and Laplace equation 
V =0. It may be noted that the expressions for J given by Smith (1962) and Tidwell 
(1970) contain errors. 

The polyhydral frame of Aryabhata with 26 faces can be closely approximated by a 
spherical shell. This approximation results in considerable simplification without 
introducing significant errors for such satellites (Tidwell 1970). For a spherical shell. 


®=constant, <o =&>k, 


fxydV—fxzdV— jyz dV—0, 

jx i dV=fy 2 dV=fz 2 dV=A7rr i sl3. (3) 



coordinotes 
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(Vernal equinox) 


(a) geometry of orbit 



Figure 1. Co-ordinate system for description of satellite motion 
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Hence, 

T=K^{B x B y i +B I B y i-(Bf+B*)k}, 

(4) 

where 

K 1 =2irr i s<x/3. 

(5) 

The component of the torque causing spin-decay is 



t,=-kmb x *+b*). 

(6) 


In terms of nodal components, satellite orientation and orbital position (figure 1) 
it is 


T^-KMBJ+B y f + BJ-{(B xn Cf +a +B yn S,f, +a )S e +B zn Cg }*]. (7) 

Idealising the earth’s field by that of a dipole placed at the centre of the earth and 
rotating with it, the nodal components can be written as (Wheeler 1967) 

B XH =K 2 {A+3(AC 2a +B S 2a )>, 

B rn =K 2 {B+3(A S ta —B C 2a )>, (8) 

B t =2 K 2 C, 


where A — S(S n , K t —~MJ2r 0 s , 

B = C € S,—S e C,C M , (9) 

C — —C t C t — fxi, 

Cl. 


Also, from Euler’s relation, the rate of change of spin is 



•o'—TJJ. 

(10) 

Hence, 

<•>'=—K<af(e, px, i, a, 8, <A), 

(11) 

where 

K = K x K 2 2 jl = <mr*s ( MJR*)*/6I . 

(12) 


It may be noted that the spin decay rate depends upon the satellite parameters, and 
its orientation and location. For inertially oriented spin-stabilised satellites, angle 4> 
varies as the orbital angle (fs angle 8 measured from the orbital normal 

undergoes long-period variation due to orbital precession, while varies due to the 
earth’s rotation and orbital precession. Thus the decay rate (equation (11)) contains 
a term with secular, long-period, diurnal and short-period variations. The mean 
daily spin-decay is given by 

<»'=—K<o(A 1 —B 1 5^+C 1 5 2 o), 


(13) 
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where ^=10-10 C £ 2 C, 2 -5 S € 2 S t \ 

Ri=9 S* —4+5 ( 2 (Stf+ 8-5) (1—1*5 S € 2 ), 


C!=S 2t St o (1—1*5 S e 2 ). 

(14) 

The decay will be minimum when 


Qa>'/d9=0 and dw7d^o=0, 

(15) 

or ^r 0 =7r/2 and 6= 0-5 tan -1 (—2 CJ—BJ, 

(16) 

while maximum decay occurs at 6 right angle to the above, 
mum rates are, respectively. 

The minimum and maxi- 

o'=-Ko [A x -5 1 /2-(5 1 2 +4C 1 2 )*/2], 

(17) 

and o'=-Ko [A t - S,/2+(5 1 2 +4C 1 2 )l/2]. 

(18) 


At all other orientations the rate will lie between these two limits. With this general 
development we proceed to consider the case of Aryabhata. 


3. Case study of ‘ Aryabhata ’ 

For Aryabhata, 

j== 51°, /?=6-97x 10 8 m (approx.), 7=100 kg m 2 . 

Equating its volume with that of a sphere, 
r=0-65 m (approximately). 

The shell is made of aluminium alloy of 1 mm thickness. The internal conducting 
structure contributing to eddy current losses can be accounted for by an appropriate 
increase in the shell thickness. However, since this increase is inversely proportional 
to the fourth power of the distance, the total change is not significant in the present 
case. The other constants in SI units are: 

<7=3-2x 10 7 (ohm-m) -1 for aluminium 

M e = 8-01 x 10 16 Weber-m, 

Therefore, from (12) 


K= 1-6737 x 10 -8 per second. 


(19) 
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It may be emphasised that one should be careful in using proper values of the para¬ 
meters depending on the units in use for evaluating K. The minimum and maximum 
decay rates (equation (17) and (18)) are, respectively, 

«/ min =-4 Kco at tp 0 =ir/2, 0=—114°, 

at ^ 0 =7r/2, 0=24°. (20) 

The bounds of spin-rate history are therefore 

w max="o exp (-0-00578 t ), 

£ °min ==a ’o exp (—0-00954 f), (21) 

where t is in days. These are plotted in figure 2, which also shows the design esti¬ 
mates and the observed values (Goel et al 1978). 

It is apparent that the present simplified analysis gives results which are closer to the 
actual observations. Knowing the orientation of the spacecraft one should be able 
to predict the actual spin rate history from (13), (14) and (19). From the comparison 
of the observed value with the predicted limits it appears that initially 0 must have 
been close to those for minimum decay. Due to orbital precession and perturbations 
they gradually change making the observed value drift away from the minimal. 
It may also be added that if one assumes the earth’s geomagnetic pole to be coincident 



Figure 2. Spin-rate history of Aryabhata 
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with the geographic north (i.e. e=0°), the analysis is considerably simplified with 
the predicted decay rates decreasing by about 5% in the present case. 


4. Conclusions 

Based on the above simplified analysis and comparison of results the following 
observations are made. 

(i) Eddy current losses are the main source of the spin decay of Aryabhata. 

(ii) A spherical conducting shell equating the outer frame is adequate for modelling 
the satellite. 

(iii) The decay is proportional to a non-dimensional parameter K, which increases 
with an increase in the size and thickness of the satellite’s conducting structure 
and decreases with an increase in the spin moment of inertia and the orbital 
radius. It is proportional to the conductivity of the shell material. 

(iv) The instantaneous rate of decay also depends upon the orbital position and 
orientation of the satellite. Components with secular, long-period, diurnal 
and short-period (orbital and twice orbital) variations are present. For life¬ 
time studies only the first two are important. Knowing the orientation and 
spin-rate at a given time it is possible to generate the entire time history. 

(v) In general, for an orbit of higher inclination the decay rate is more. 

It is felt that the above analysis is applicable to many spin-stabilised satellites. 


Thanks are due to Dr A S Prakasarao and Mr P S Goel of the ISRO Satellite Centre 
for providing the satellite data used in this paper and for useful discussions. 


List of symbols 


B 

c v s q 

I 


i,j,k 

J 

K, K lt K 2 

M. 

R 


r 

s 

T 

t 


magnetic field intensity 
cos ( q ), sin (q) 

satellite’s moment of inertia about spin axis 

orbital inclination 

unit vectors along x, y, z 

eddy current intensity 

constants 

earth’s magnetic dipole moment 

orbital radius 

satellite’s equivalent radius 

satellite’s shell thickness 

torque 

time in days 
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V volume of conducting material 

x, y, z principal co-ordinates 

a position in orbit, measured from ascending node 

e angle between the earth’s geographic and geomagnetic north 

p position of line of intersection of geomagnetic and geographic equator 

measured from Vernal equinox (figure 1) 
to, u> satellite spin rate, spin decay rate, respectively 

£2 argument of ascending node of the orbit (figure 1) 

distance from the centre of mass (S) to an elemental volume dV 
ip, 0, <p Euler’s angle defining principal axes S-xyz with respect to orbital frame 
S-x 0 y a z 0 (figure 1) 

Subscripts jc, y, z and xn, yn, zn denote components along principal and nodal co¬ 
ordinates, respectively. 
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